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Direct detection of Milky Way wimps are discussed within the framework 
of R-parity conserving SUGRA models with grand unihcation at Mq- Two 
questions are discussed: what SUGRA models can account for the DAMA 
data if this data is conhrmed, and is the full SUGRA parameter space ac¬ 
cessible to future planned detectors. Experimental constraints of the Higgs 
mass bound, the 6 —>• sy bound, relic density constraints (including all co¬ 
annihilation channels), etc. are imposed. In addition, the effect of the pos¬ 
sible muon g — 2 anomaly are examined. For mSUGRA, we hnd that the 
Higgs mass and b ^ sj constraint puts a lower bound mi /2 > (300 - 400) 
GeV (i.e. m^o > (120 - 160) GeV) for tan/? <50, and thus the largest theo¬ 
retical neutralino-proton cross sections still lie signihcantly below the DAMA 
3a lower bound. (Predictions for tan/? > 50 become sensitive to the precise 
value of mt and m^.) If in addition one imposes the muon anomaly constraint, 
/i must be positive and an upper bound of mi /2 < 850 GeV for tan/? < 50 
is obtained. More generally, if /i > 0 and mi /2 < 1 TeV, the cross sections 
are ~ 10“^*^ pb, and hence this parameter space would be mostly accessible 
to planned high sensitivity detectors. For non-minimal SUGRA models, the 
cross sections can be considerably larger, and a simple SU(5) model with non- 
universal soft breaking in the Higgs and third generation is seen to give cross 
sections in the DAMA range for tan /? ~ 15 with m^o > 80 GeV, and minimum 

cross sections ~ 10“^° pb for /r > 0. 
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1 Introduction 


It is generally agreed that a large part of the Milky Way galaxy is made of dark matter, 
i.e. current estimates of the dark matter energy density are pdm — (0.3 — 0.5)GeV/cm^. 
However, what this dark matter is made of is much less clear, and many suggestions 
exist. We will here assume a SUSY explanation, and in particular assume that the dark 
matter particle is the lightest supersymmetric particle (the LSP), the lightest neutralino, 
Xi- A number of ways have been suggested to detect SUSY wimps: One may look for 
annihilation in the halo of the Galaxy leading to e^, p + X and search for 

the anti-particle signals. Alternately, one expects neutralinos to be captured by the Sun 
or Earth, fall into their center and there annihilate. Muon neutrinos could then escape, 
and one could look for the signal + Xj —>■ X with neutrino telescopes on Earth. 
Finally, there is the possibility of direct detection by the scattering of incident neutralinos 
by terrestrial nuclear targets. 

The theoretical expectations for these possible signals depends upon the model chosen. 
The most general SUSY model, the MSSM, contains over 100 free parameters (63 real 
parameters!), and so it does not have great predictive power. The supergravity (SUGRA) 
models with grand unihcation at the GUT scale Mq ~ 2 x 10^® GeV Bi, apply to a 
wide range of phomenena, and are relatively predictive. We will consider here then the 
SUGRA models, both those with universal soft breaking (mSUGRA) and non-universal 
soft breaking at the GUT scale. Our models will also possess R-parity invariance, to 
guarantee the existence of a stable LSP dark matter candidate. In particular we examine 
the gravity mediated models. (Other possibilities such as gauge mediated soft breaking 
have difficulty in constructing a viable dark matter candidate, while the anomaly mediated 
models appear to be in contradiction with the recent Brookhaven muon magnetic moment 
data.) 

Goncerning possible signals, there has been indications of an excess of p and events 
in the halo, indicating possible halo x\ annihilations. This data is quite interesting and has 
led to some recent analyses P, H, though it is difficult to be sure of all the astrophysical 
backgrounds. Analyses have also been given of the possibilities for present or future 
neutrino telescopes (AMANDA, Ice Gube, ANTARES) to see energetic from the Sun 
1^, ^ . Thus 1^ hnds that Ice Gube and ANTARES can be sensitive to Xr annihilation for 
non-universal SUGRA models but only for large tan (3 i.e. tan [3 ~ 35. The most promising 
way of detecting halo wimps for a wide range of parameters remains, then, direct detection 
by terrestrial targets, and we will restrict the discussion to this possibility from now on. 

Experiments for direct detection of Milky Way wimps is now entering a new (and 
exciting) phase. There exists the DAMA results with data that indicates the observation 
of the annual modulation signal (at the 4 -ct level), and the GDMS exclusion curves which 
appear to exclude much of the allowed DAMA region. However, there are a number of 
new experiments that should be able to clarify matters in the relatively near future. Thus 
DAMA is constructing an upgrade to a 250 kg Nal detector (from their current 95 kg 
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detector). Two other detectors that will be completed in the near future, GENIUS-TF 
and ZEPLIN II, should be able to give independent observations of whether the annual 
modulation effect exists, and when CDMS is moved to the SOUDAN mine, its exclusion 
curves will be sensitive to the full DAMA allowed region. Further, in the more distant 
future, a number of new detectors are being planned with much higher sensitivity, e.g. 
GENIUS, Gryoarray, ZEPLIN IV, and GUORE. 

In view of this, we consider here two questions: (1) If the DAMA data is conhrmed, 
what SUGRA models can account for the relatively large neutralino-proton cross sections 
this data implies, and (2) what are the lowest SUGRA cross sections theoretically pre¬ 
dicted, and will the planned detectors be able to scan the full SUGRA parameter space. 
We will see that the DAMA data, if conhrmed, would exclude the mSUGRA model (but 
be consistent with non-universal models), and if the Brookhaven g^ — 2 anomaly is valid, 
future detectors will be able to test SUGRA models over the full parameter space and 
hence be competitive with accelerator experiments. 


2 Experimental and Theoretical Constraints 

One of the important features of SUGRA models is that they apply to a wide range of 
phenomena, and data from these already restrict the parameter space and thus strengthen 
the predictions. We list here the main experimental constraints that we use. 

(1) Higgs mass 

The LEP data 0 places the lower bound of rrih >114 GeV. However, theoretical cal¬ 
culation of the light Higgs mass, rrih, still have an uncertainty of about 3 GeV. Thus 
we will conservatively interpret the experimental bound to require that the theoretically 
calculated mass obey rrih >111 GeV. 

(2) 6 ^ sy branching ratio 

The GLEO data has both systematic and theoretical uncertainties. We therefore take 
a relative broad range around the GLEO central value of 

1.8 X 10“^ < B{B X, 7 ) < 4.5 x 10"^ (1) 


(3) Xi relic density 

The cosmic microwave background radiation (GMB) and other data now give a rela¬ 
tively constrained value for the cosmic mass density of dark matter. This is measured in 
terms of the value of Here Ddm = p/ Pc where p is mass density of dark matter, 

and pc = /SirGjsf is the critical density to close the universe {H = [lOOkm/sMpc)h 
is the current Hubble constant. Gat is the Newton constant). A recent analysis gives 
= 0.139 ± 0.026 p. We take a 2.5cr range around the mean: 

0.07 < < 0.021 (2) 
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In looking for the minimum cross section, we conservatively extend this range to be from 
0.025 to 0.25. 

(4) Recently, the BNL E821 experiment has observed a 2.6cr anomaly in the muon 
magnetic moment — 2)/2. Such an anomaly was predicted some time ago for 

mSUGRA models O, O and we will thus interpret the affect as being due to SUGRA. 


We take here a 2a range around the E821 central value 


11 X 10"^° < < 75 X 10“^° 

— — 


(3) 


While this effect is not yet certain (and there is some debate as to the size of the Standard 
Model hadronic contribution, see e.g. |]^), it has a signihcant impact on dark matter 
prediction. We will therefore include it in our analysis, but show explicitly what effects 
it produces. Further data, currently being analyzed by E821, should reduce both the 
experimental error and theoretical uncertainty by more than a factor of two, and results 
should be available early next year Q 

In order to carry out accurate calculations, it is necessary to include a number of 
theoretical corrections, and we list some of these here. (1) Two loop gauge and one loop 
Yukawa renormalization group equations (RGE) are used from Mq to the electroweak scale 
Mew and QGD RGE are used below for the light quark contributions. We chose Mew = 
(2) Two loop and pole mass corrections are included in the calculation of 
ruh- (3) One loop corrections to rrih and [[ 
large tan/?. (4) Large tan/? NLO SUSY corrections to b 


15[| are included which are important for 
pT| , p!7| are included. (5) 
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All stau-neutralino co-annihilation channels are included in the relic density calculation, 
this analysis being done for both small and large tan/? |^, We note that we do 

not include Yukawa unihcation or proton decay constraints, as these depend sensitively 
on post-GUT physics about which little is known. (For example, string theory with e. g. 
an SU(5) GUT group requires grand unihcation of the gauge coupling constants, just as 
SUGRA theory does, but the Yukawa or proton decay constraints of SUGRA need not 
apply ig.) 


3 mSUGRA Model 


The minimal supergravity model, mSUGRA, is the most predictive of the models as it 
depends on only four additional parameters and one sign. It is convenient to chose these 
parameters as follows: mo, the scalar soft breaking mass at Mg] mi/ 2 , the gaugino soft 
breaking mass as Mg- (We note here the approximate formulae that m^o = 0.4mi/2, 


*Recently, two papers have appeared (Knecht et.al. ( hep-ph/0111059| ), Hayakawa and Kinoshita ( |hep- 
ph/0112102| )) implying a sign error in the scattering of light by light contribution to a^. This reduces 
the 2.6(7 effect to 1.6a, our lower bound then becomes la below the central value. We note that with the 
further data, the E821 experiment would be expected to see an ~ 4cr effect, if the reduced central value 
were to remain unchanged 
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the chargino mass m-± = 0 . 8 mi/ 2 , and the gluino mass rrig = 2.5mi/2.); the cubic 
soft breaking mass at tan/? =< H 2 > / < Hi > a.t the electroweak scale (where 
hf( 2 ,i) gives rise to up,down quark masses); and the sign of /r, the Higgs mixing parameter 
(which appears in the superpotential as fiHiH 2 ). We examine these parameters over the 
following range: mo, mi /2 < 1 TeV, 2 < tan/? < 50, and Aq < 4 mi/ 2 . The range of mg 
and mi/2 is sufficient to explore the full range that will be accessible to the LHC (and we 
will see that if the — 2 anomaly is real, it will cover the full parameter space). There 
is a small additional allowed region of very high tan /? i. e. tan /? = 50 — 55. However, as 
we will see below that predictions in this region are very sensitive to the precise values of 
the t and b quark masses which are not at present known. 

In general, the neutralino-nucleus scattering cross section has a spin dependent part 
and a spin independent part. However, the transition from quark to proton scattering 
for heavy nuclei, the spin independent cross section dominates, and this allows one to 
extract (to a good approximation) cr^o_p, the neutralino-proton cross section, from any 
dark matter detector data. The basic quark diagrams for neutralino-quark scattering, 
crj^o_q, are s-channel scattering through hrst generation squarks, and t-channel scattering 
through h and H (where H is the heavy neutral Higgs boson). One can determine cr^o_p 
then from the basic quark scattering. However, this can involve errors perhaps of a factor 
of ~ 2 from uncertainties in the strange quark content of the proton and the pion-nucleon 
cr term, CTr-N- (We use here the most recent evaluations of the a term, (7.,t-N = 65 MeV 
Earlier evaluations, a-n-N = 45 MeV,would reduce the cross sections by perhaps 
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a factor of 2 .) 

The parameter space for calculations of the above detector cross section are con¬ 
strained signihcantly by the requirement that the theory predict the right amount of relic 
neutralino dark matter. In the early universe, two neutralinos can annihilate to Standard 
Model particles through s-channel Z, h, H and A poles {A is the CP odd Higgs boson), 
and through t-channel sfermion diagrams. However, if a second particle becomes nearly 
degenerate with the Xi^ must include it in the early universe annihilation processes, 
leading to the phenomena of co-annihilation. In SUGRA models with universal gaugino 
masses at the GUT scale, this “accidental” near degeneracy occurs naturally for the light 
stau, fi. To see this semi-quantitatively, we note that for low and intermediate tan/?, the 
RGE can be solved for the right selectron, cr, and x^ masses to give at the electroweak 
scale: 


rriej^ = 1^11 + 0.15m^/2 ~ sin^OwM"^ cos 2/? (4) 

m^o = 0.16m^/2 

where the numerical factors come from solving the RGE. The last term is about (35GeV)^, 
and so for mo = 0, the cr becomes degenerate with the at mi /2 = 350 GeV, and thus 
co-annihilation begins at mi /2 = (350 — 400)GeV. As mi /2 increases, mo must be raised 
in lock step (to keep mg^ > m^o). More, precisely, it is the light stau, which is the 
lightest slepton, that dominates the co-annihilation phenomena. However, one ends up 
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Figure 1: Corridors in the mo — mi /2 plane allowed by the relic density constraints for 
(bottom to top) tan/3 = 10, 30, 40, Aq = 0 and /i > 0. The lower bound on mi /2 is due to 
the ruh lower bound for tan/3 = 10, due to the b ^ s'f bound for tan/3 = 40, while both 
these contribute equally for tan/3 = 30. The short lines cutting the channels represent 
upper bound from the 5 '^ — 2 experiment. [24] 


with corridors in the mo —mi /2 plane with allowed relic density, with mo closely correlated 
with mi /2 and mo increasing as mi /2 increases. 

We now turn to the other experimental constraints. The Higgs mass bound and the 
b ^ sj bounds produce comparable constraints on the parameter space, the Higgs mass 
generally being dominant for low tan /3 and the 6 —> sy for high tan /3. Over the entire 
tan/3 range these constraints produce a lower bound on mi/ 2 : 

mi /2 > (300 - 400)GeV (5) 

corresponding to a lower bound on the neutralino mass of 

m^o > (120 - 160)GeV (6) 

This means that most of the parameter space is in the f — Xi co-annihilation domain 
in the relic density calculation. Thus relic density bounds then imply ttiq is approximately 
determined by mi/ 2 . 

If we now include the Brookhaven — 2 constraint we obtain additional restrictions 
on the parameter space. First the sign of implies that /i > 0. Further, the lower 

bound on implies an upper bound on mi/ 2 . Thus for Aq = 0, one has mi /2 < 440 

GeV for tan/3 = 10, mi /2 < 790 for tan/3 = 40 and mi /2 < 850 for tan/3 = 50. Thus 
assuming the Brookhaven E821 data, we obtain both an upper bound on mi /2 and a lower 
bound on mi/ 2 - The reduction of the parameter space for small tan/3 then puts a lower 
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Figure 2: Corridors in the mo — mi /2 plane allowed by the relic density constraint for 
tan/3 = 40, /i > 0 and (bottom to top) Aq = 0, —2mij2, 4mi/2. the curves terminate at 
the lower end due to the b ^ s'y constraint except for/lo = 4mi/2 which terminates due 
to the ruh constraint. The short lines cutting the corridors represent the upper bound on 
mi /2 due to the — 2 experiment. [25] 


bound on tan [3. We hnd 


tan/3 > 7(5), for Aq = 0(—4mi/2) 


(7) 


Some of the above results are illustrated in Figs. 1 and 2. Fig. 1 shows the allowed 
regions (which are the co-annihilation channels) in the mo — mi /2 plane for tan/3 = 10, 
30, 40 (from bottom to top) for Aq = 0 The lower mi /2 bounds are due to the rrih 
bound for tan/3 = 10, and due to the b ^ S'j bound for tan/3 = 40. (They are equally 
constraining for tan/3 = 30.) The g^ — 2 upper bounds are given by the short lines cutting 
the channels. One sees that the co-annihilation corridors occur for higher mo as tan/3 
increases. Fig. 2 shows the Aq dependence for the case of tan/3 = 40. One sees that the 
co-annihilation corridors he higher in mo for Aq different from zero. 

For tan/3 > 50, a new phenomena arises due to the fact that the heavy Higgs (H, A) 
become light, and can contribute to the annihilation cross section. Results in this domain 
become very sensitive to the precise values of rrit and mb @]. We illustrate this in Figs. 
3 and 4. Fig. 3 gives the allowed region in the mo — mi /2 plane for tan/3 = 50, Ao = 0 
using the central experimental values of m^ = 175 GeV and m^ = 4.25. However, if mt is 
reduced by la to 170 GeV, one sees from Fig. 4, that the shape of the allowed region is 
signihcantly modihed, and a new corridor is opened at low mi /2 rising to mo > 1 TeV. 

We see that the parameter space has become highly restricted, both at the lower 
end and the upper end, and this should help in clarifying the predictions for dark matter 
detection. The general features of a^o_p that explain its properties are that a^o_p increases 
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Figure 3: Allowed region in the mo — mi /2 plane for tan/3 = 50, Aq = 0, /i > 0 for 
rrit = 175 GeV, mb = 4.25. The lower and upper bounds on mi /2 are due to the b —>■ s'} 
and ttfj, constraints respectively. 


with increasing tan/3, and decreases with increasing mi /2 and mo (and as we’ve seen,mi /2 
and mo move together). Thus the maximum values of cr^o_p are generally expected to 
occur for high tan/3, and low mi/ 2 , mo. We now return to our two questions of Sec. 1: 
in mSUGRA, how large can cross sections get, i. e. can they accommodate the DAMA 
data, and how small can they become, i. e. will future detectors be able to scan the entire 
parameter space. The current DAMA annual modulation signal is in the range 
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(o'wimp)DAMA = (10 - 10 )pb and rn^imp < (100 - 300)GeV 


( 8 ) 


where the uncertainty in the upper bound on the wimp mass being due in part to as¬ 
tronomical uncertainties about the Milky Way. We compare this now to the mSUGRA 
predictions. Fig. 5 shows the Aq dependence of the cross sections for tan/3 = 40. Note 
that Ao = 0 gives the largest cross section. We see that it is the low mi /2 bound that 
prevents the cross section from becoming large, (and in fact if the bound on mh were to 
increase to 120 GeV, the lower bound on mi /2 would rise to 200 GeV for Aq = —2 mi/2, 
to 215 GeV for Aq = 0, and to 246 GeV for Aq = 4mi/2 further reducing the maximum 
cross sections). Fig. 6 shows the cross sections for tan/3 = 10, 30, 40 and 50 for Aq = 0, 
mt = 175 GeV. The largest cross sections occur for the largest tan/3 and smallest mi/ 2 , 
as expected. We hnd in general, 

(^x\-v ^ (0-07) X 10“®p6 for tan/3 < 50 (9) 

which is well below the DAMA bound. Thus it is not possible for the mSUGRA model to 
accommodate the DAMA data. We note that central in this result is the mh and b ^ s'^ 
constraints which put a lower bound on mi /2 and hence require m^o to be relatively high. 
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Figure 4: Allowed regions in the mo — mi /2 plane for tan/3 = 50, Aq = 0, /i > 0 for 
rrit = 170 GeV, mb = 4.25 GeV. The lower and upper bounds on mi /2 are due to the 
b ^ sj and constraints respectively. 


We turn now to the second question, what are the minimum values possible for o'^o-p, 
and will future detectors being planned be able to observe them. Minimum cross sec¬ 
tions arise from low tan/3 and large mi/ 2 , mo. Here the question of the validity of the 
Brookhaven ( 7 ^ — 2 data plays an important role in two ways. First this data would elim¬ 
inate the p < 0 possibility and for p < 0, there are accidental cancellations that can 
reduce cT^o_p to well below 10“^° pb over wide ranges of parameters [|^ |2^, which would 
make this part of the parameter space experimentally unaccessible. Thus these extremely 
low cross sections would be eliminated. Second, it is the lower bound on that 

produces the upper bound on mi/ 2 , and this effect is larger for smaller tan/3 (see e. g. 
Fig. 6) moderating the reduction of the cross section as tan/3 decreases. Fig. 7 shows 
(T^o_p for tan/3 = 10, Aq = 0 (upper curve), —4mi/2 (lower curve). We see here that 
a. 


X/ 


o_p > 4 X 10 pb. The parameter space ends for tan/3 = 7 (5) for Aq = 0(—4 mi/2) 


and in general one hnds cr^o_p >10 pb. 

Future detectors hope to obtain a sensitivity of o'^o_p = (10“® — 10“^°) pb. Thus 
they should be able to scan most of the mSUGRA parameter space. If this sensitivity is 
achievable, dark matter experiments would be competitive with accelerator experiments. 


4 Non-Universal Models 

We next investigate the possibility that the DAMA data gets confirmed, and the large 
neutralino -proton cross sections it implies are correct. To accommodate this theoretically, 
one can consider models with non-universal soft breaking at Mq in the Higgs and third 
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Figure 5: c^^o-p for tan/? = 40, /i > 0, for (bottom to top), Aq = — 2 mi/ 2 , = 4 mi/ 2 , 

Aq = 0. The lower bounds on m^o are as in Fig. 2. 


generation squarks and sleptons. A convenient parameterization is 


rriul 


= ml{l + 6i); = ml{l + 62); 


= rrin 


(1 + ^ 3 ); m^l = ml{l + Si); m = ml{l + S 5 ); 


= mo(l + ^ 6 ); rriif = ml{l +67). 


( 10 ) 


with — 1 < < 1. While this appears to involve introducing a large number of new 

parameters, one can understand the effects they produce relatively simply, as much of 
the physics of what occurs is governed by the fi parameter. One can qualitatively see the 
effects of the non-universalities for low and intermediate tan/? where one may solve the 
RGE for /i^ analytically PP| : 


/i^ = 




- 1 
1 + -Do 


A-3Do , 1, , 1-Do,, , ,, 


A 


tUq-|- universal parts -|- loop corrections. 


( 11 ) 


where t = tan/? and Ho = 1 — (^t/ 200 sin2/?)2. One see that Ho is small (i. e. Ho = 0.25) 
and that the universal contribution to the term nearly cancels out. Thus one does 
not need a large amount of non-universality for the Si contributions to dominate. Now 
the important point is that if decreases (increasing the Higgsino content of y?) o'^o_p 
increases, and similarly if increases, o'^o_p decreases (though by not as large an amount). 
From Eq(ll), then we see that the cross section increases if ^ 1 , 3,4 < 0, ^2 > 0. We consider 
here a simple SU{5) model where (5io = 53 = 54 = 55 , and S 5 = Sq = Sj. We chose 
5i = —1, S 2 = 1, with 55 > 0 and 5io < 0. With these choices, is reduced and hence 
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Figure 6: o'^o_p for Aq = 0, m* = 175 GeV, rrih = 4.25, for tan/? = 10, 30, 40 50 (from 
bottom to top). The lower bound on m^o is due to the h ^ sj constraint for tan/3 = 30, 
40, 50, and from the constraint for tan/3 =10. 



the cross section is increased. One must, however, make sure the relic density constraint 
is simultaneously satished. This occurs due to the fact that the <5, chosen also lowers uia- 
Thus the RGE for rriA give for low and intermediate tan/3: 


ml 


r + 1 
l + Dq 


3(1 


Do) I 1 


D„ 


52 + (3i 


(^3 + 64 


tUq + universal parts + loop corrections. 


( 12 ) 


Thus mA is simultaneously lowered opening an annihilation channel for through 
an A s-channel pole, and allowing the relic density constraint to be satished. What is 
happening is similar to what happens in mSUGRA at tan /3 > 50, but for the non-universal 
model this effect can happen at low tan/3 and with a signihcantly larger effect. Fig. 8 
shows the maximum value of o'^o_p for the model discussed above for tan/3 = 12 (lower 
curve), 18 (upper curve). One sees that one can generate o'^Q_p in the DAMA domain 

for tan/3 ~ 15. Note also that the lower bound on m^o is greatly reduced, since the 
Higgs mass and b ^ constraints are much less conhning for the non-universal model. 
(Again a lower wimp mass would be in better accord with the DAMA results.) Thus if the 
DAMA data is conhrmed by other groups, this could point to the validity of non-universal 
SUGRA models. 

The minimal cross sections arise from low tan/3, large mi/ 2 , and also from reversing 
the signs of the 5i. One hnds for /r > 0 (as implied by the g^—2 data) that cr^o_p > 10“^°pb, 
which as before is within the reach of future planned detectors. 
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Figure 7: o'xi-p tan/3 = 10, Aq = 0 (upper curve), Aq = —4mi/2 (lower curve). The 
lower bound on m^o is due to the mt bound for Aq = 0, and the b ^ S'j constraint for 
Ao = -4mi/2. 

5 Conclusions 

We have discussed here the direct detection of neutralino dark matter in the Milky Way 
within the framework of R-parity conserving gravity mediated supergravity models with 
grand unihcation at Mq = 2 x 10^® GeV. In particular, we examined two questions: if 
the DAMA data is confirmed, can SUGRA models accommodate the large cross sections 
implied by the data, and also, what are the lowest cross section the SUGRA theory 
predicts and will future planned detectors be able to scan the entire parameter space. 

There are now many experimental constraints that allow one to make clearer predic¬ 
tions for these models. Thus for mSUGRA, the Higgs mass and 6 —>■ sy constraints create 
a lower bound of mi /2 > (300 — 400) GeV (i.e. m^o > (120 — 160)GeV ), which puts the 
parameter space mainly in the co-annihilation domain, and (approximately) determines 
mo in terms of mi /2 (for hxed tan/3 and Aq)- We hnd then that the largest cross section 
is signihcantly below the lower bound of the DAMA data for tan (3 < 50. (For tan (3 > 50, 
predictions are sensitive to the precise values of and mi,.) If the muon g — 2 anomaly 
is conhrmed by future data, then significant further reduction of the parameter space 
occurs. One has that p is positive, eliminating the possibility of very small cross sections. 
In addition, this effect produces an upper bound of mi /2 < 850 GeV for tan/3 < 50 as 
well as a lower bound on tan/3 of tan/3 >7 (5) for Ao = 0(—4mi/2)- We find then over 
the remaining parameter space that cross sections are greater than ~ 10“^° pb and hence 
most of the parameter space will be accessible to the future planned experiments. This 
result also holds without the — 2 anomaly for the /i > 0, mi /2 < 1 TeV part of the 
parameter space. 

For non-minimal SUGRA models, with non-universal soft breaking in the Higgs and 
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Figure 8: Maximum values of cr^o_p for tan/? =12 (lower curve), 18(upper curve) for 

^2 = 1 = ^10 < 0 ) > 0 ) ^0 = 0 , /i > 0 . 

third generation squark and slepton sectors, one can generate models with significantly 
larger neutralino-proton cross sections. Thus for a simple SU(5) model, one finds it is 
possible to find cross sections in the DAMA region for tan/? > 15, and also the lower 
bound on the neutralino mass can be reduced to 80 GeV. These models satisfy the relic 
density constraint by simultaneously reducing the A mass (even for the relatively low 
tan/? above), and one also finds that for the /r > 0 part of the parameter space the 
cross sections should not fall significantly below 10“^° pb, and thus be experimentally 
accessible. 
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